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ABOUT THE EMISSION LINE SEQUENCE OF H II GALAXIES 

Grazyna Stasinska 1 , Yuri Izotov 2 
RESUMEN 

El resumen sera traducido al espanol por los editores. 

We consider 400 H n galaxies in which the oxygen abundances were obtained by electron temperature 
based methods. We split the sample in three metallicity bins. In each bin, a narrow sequence is found not only 
in pure emission line ratio diagrams but also in terms of H/3 equivalent width, indicating the existence of an 
evolutionary sequence. Our diagrams show unambiguously the existence of an evolution on a timescale of a few 
Myr. We compare the observed sequences with photoionization models of appropriate metallicity. In order to 
understand the evolution of H n galaxies one needs to consider the evolution of the gas as well as that of the 
stars. The observations require EW(H/3) to decrease more rapidly than predicted by the passive evolution of a 
star burst. A photoionized adiabatic expanding bubble with a covering factor decreasing with time reproduces 
most diagrams. However, the question of the heating of H n galaxies and the origin of the nebular He n A4686 
emission are not settled. We find evidence for self-enrichment in nitrogen on a time scale of several Myr. 

ABSTRACT 

We consider 400 H n galaxies in which the oxygen abundances were obtained by electron temperature 
based methods. We split the sample in three metallicity bins. In each bin, a narrow sequence is found not only 
in pure emission line ratio diagrams but also in terms of H/3 equivalent width, indicating the existence of an 
evolutionary sequence. Our diagrams show unambiguously the existence of an evolution on a timescale of a few 
Myr. We compare the observed sequences with photoionization models of appropriate metallicity In order to 
understand the evolution of H n galaxies one needs to consider the evolution of the gas as well as that of the 
stars. The observations require EW(H/3) to decrease more rapidly than predicted by the passive evolution of a 
star burst. A photoionized adiabatic expanding bubble with a covering factor decreasing with time reproduces 
most diagrams. However, the question of the heating of H n galaxies and the origin of the nebular He II A4686 
emission are not settled. We find evidence for self-enrichment in nitrogen on a time scale of several Myr. 
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1. INTRODUCTION 

It has been known for years that giant H n regions 
form a narrow sequence in various emission line ratio 
diagrams (e.g. McCaU, Rybski & Shields 1985). The 
reasons for the existence of such a sequence are not 
yet fully understood. The first interpretations, based 
on single star photoionization models, concluded that 
metallicity is driving the sequence, but variation of 
an additional parameter - either the effective tem- 
perature of the stars or the ionization parameter - 
is needed to reproduce the observed sequence (Mc- 
Call ct al. 1985; Dopita & Evans 1986). More recent 
studies using stellar population synthesis have shown 
that the spectral energy distribution of the ioniz- 
ing radiation depends naturally on metallicity, both 
through stellar structure and evolution and through 
the effect of opacities in the stellar atmospheres (e.g. 
Mc Gaugh 1991). 
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However, star formation proceeds by bursts and 
giant H II regions are powered by clusters of co- 
eval stars (e.g. Sargent & Searle 1970, Mas-Hesse & 
Kunth 1999). Stellar evolution produces a gradual 
modification of the integrated stellar energy distri- 
bution. Spectra of H II regions whose ionization is 
dominated by one single cluster can then be caught 
at different cluster ages. This is the case of blue 
compact dwarf galaxies or dwarf irregular galaxies, 
which are referred to under the common name of H II 
galaxies. These objects make up invaluable tools to 
study the evolution of stellar systems and their in- 
teraction with the surrounding interstellar medium 
during the phase where massive O stars are present. 

Several studies (see references in Stasinska & Izo- 
tov 2003, hereafter SI2003) have been devoted to 
this question, based on various samples and some- 
what different approaches. Here, we give an account 
of the SI2003 paper, which aimed at finding what 
conditions are needed to reproduce the observed se- 



1 



2 



STASINSKA & IZOTOV 




10 100 
EW(H0) (A) 



412 



-0.1 



[Oil] 



10 100 
3727 / H/S 

4o'l" 



10 100 
EW(H0) (A) 



: b 


412 ; 












...I I 



1000 



10 

EW(H0) (A) 



e 


392 ; 







^O -3 0.01 0.1 1 
[Ol] 6300 / H/3 



o 




. , 48 ; 


X ^_ 

\ 9 
o 

CD 
CO 

to 






= o 

CD *~~ 
X 











i 



10 100 
EW(H0) (A) 



c 


392 \ 












.1 



o 
+ 



10 100 
EW(H/9) (A) 



1000 



: f 




"I ! 

412 ; 
















Fig. 1. The entire sample of H II galaxies in diagrams involving line ratios and EW (H/3). The number of data points is 
shown in each panel. Byurakan objects are in red, SDSS objects are in blue. 



quence of H n galaxies with models of H II regions 
surrounding evolving starbursts. The success of this 
enterprise has greatly benefited from the fact that 
we were able to construct a homogeneous sample of 
H II galaxies of unprecedented size, in which the oxy- 
gen abundances were derived in a model-independent 
way. This allowed us to divide the sample into three 
metallicity bins, and study each of them indepen- 
dently. 

We describe the data sample, comment on the 
observational trends and present some model se- 
quences. Details about the sample and the mod- 
cling procedure can be found in SI2003. We take the 
opportunity of this conference to develop on some 
aspects that were only briefly mentioned in SI2003. 
In particular, wc refer to the viewpoint presented by 
R. Terlevich at this conference. 

2. THE OBSERVATIONAL SAMPLE 

Ideally, the sample should be as much as possible 
relevant, homogeneous, complete (or with well un- 



derstood biases) and large (since there are at least 
two independent parameters: age and metallicity). 

Metal-poor H II galaxies are the best candidates 
to study the evolution of giant H II regions. As op- 
posed to giant H II regions in spiral galaxies, their 
spectra are not affected by the light from the stars 
located in the disk. Also, the high enough electron 
temperature allows detection of the [O in] A4363 
line and therefore direct determination of the oxy- 
gen abundance by the electron temperature method 
(this is not the case of giant H II regions found in the 
inner parts of disk galaxies). In addition, the dust 
content in H II galaxies is small, which makes the 
interpretation of the data easier and more robust. 

Our sample of H II galaxies was obtained by 
merging ~ 100 blue compact galaxies from the First 
and Second Byurakan surveys and ~ 300 emission- 
line galaxies from the early data release of the 
SLOAN digital sky survey (SDSS, Stoughton et al. 
2002). Only objects with a detected [O in] A4363 
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Fig. 2. The "high" metallicity bin. The data points correspond to the objects of our sample for which the estimated 
error in the [O in] A4363 intensity is less than 50% and which have O/H > 10~ 4 . The number of data points present 
in each diagram is shown in the upper right of each panel. Overplotted are evolutionary sequences of models with 
metallicity Z — 0.2 Z©and constant density. Each sequence corresponds to a different initial ionization parameter. 



emission line were included. The spectra of the en- 
tire sample were reduced and reddening-corrected 
with identical procedures and the abundances were 
derived in the same way. 

The spectra of the Byurakan survey galaxies 
were obtained by Izotov and coworkers to deter- 
mine the pregalactic helium abundance, therefore 
the signal-to-noise is high and the H/3 equivalent 
widths, EW(R0), arc large by observer's selection. 
The SDSS galaxies reach larger distances, but their 
H/3 luminosities are rougly similar to those of the 
Byurakan galaxies. For the majority of the objects, 
the extinction is small (E(B - V) < 0.2). 

3. DIAGNOSTIC DIAGRAMS 

The diagrams were carefully chosen to provide 
strong constraints on the evolution of the ionizing 
clusters and of their surrounding H II regions. The 
relatively high quality of the data for the entire sam- 



ple allowed us to use not only the classical strong 
lines, i.e. [O n] A3727, [O III] A5007, and H/3, but 
also weaker lines such as [O i] A6300, He I A5876 and 
He II A4686 (the latter, being at the level of 0.003 - 
0.03 of H/3, is detected only in the Byurakan sample). 
These weak lines should not be discarded as they 
provide very important diagnostics. Wc also con- 
sider a diagram involving [N n] A6584/[0 n] A3727. 

Fig. 1 shows the distribution of the observational 
points in 9 diagrams. The relative distribution of the 
Byurakan objects (red points) and of the SDSS ob- 
jects (blue points) implies that the nature of the ob- 
jects is likely similar and justifies the merging of the 
two samples. Fig. 1 shows that H II galaxies form a 
sequence not only in the classical line-ratio diagrams 
([O in] A5007/H/3 vs [O n] A3727/H/3 and [O in] 
A5007/H/3 vs [O i] A6300/H/3 3 ) but also, and very 



the sequence in the pure line-ratio diagrams is truncated 
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Fig. 3. The "high" metallicity bin. Same data points as in Fig. 2. Overplotted are evolutionary sequences of expanding 
bubble models with metallicity Z = 0.2 Zq. The green curve represents the pure expanding bubble sequence. The blue 
curve is obtained by adding an underlying older population. The red curve is obtained from the blue one assuming a 
covering factor decreasing with time. 



clearly, as a function of EW(H(3). This indicates an 
evolutionary sequence: either on a timescale of a few 
Myr, as the death of massive stars progressively re- 
duces the H/3 emission or on time scales of Gyr, as 
previous generations of stars gradually build up the 
H/3 continuum. The latter view is the one adopted 
by Terlevich in these proceedings, using different ar- 
guments 4 . The fact that the [O ill] A5007/[O n] 
A3727 and [O i] A6300/QO n] A3727+[0 m] A5007) 
ratios show very strong trends with EW(Hf3) and 
that the same trends are seen for each metallicity 
bin (see Figs. 2 and 4 and SI2003 for more details) 
argue for the first interpretation as being the main 
cause of the observed sequence in our sample. There- 



with respect to diagrams built with H II regions from spiral 
galaxies, since the metal-rich objects are lacking. 

4 The sample used by Terlevich contains also metal-rich ob- 
jects where the effects of previous generations of stars are ex- 
pected to be stronger. 



fore, EW(Hf3) is indicative of the age on the ionizing 
stellar population (although not directly, since sev- 
eral factors must be taken into account, see later). 

The fact that a sample of 400 objects shows such 
clear trends is remarkable. This means that, in spite 
of the complexity revealed by detailed studies of a 
few H II galaxies, the overall evolution of such objects 
must obey simple laws. 

Note that [N II] A6584/[0 n] A3727 is seen to 
increase as EW(R(3), although with higcr dispersion 
than [O i] A6300/QO n] A3727+[0 in] A5007). 

4. THE MODELING STRATEGY 

All the details arc published in SI2003. The ob- 
servational sample is subdivided in 3 metallicity bins: 
"high" metallicity (between 0.25 and 0.12 times so- 
lar), "intermediate" (between 0.12 and 0.036 times 
solar) and "low" (below 0.036 times solar). For each 
bin we construct sequences of photoionization mod- 
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els that represent the evolution of a giant H II region 
of appropriate metallicity. The aim is to find what 
conditions are needed to reproduce simultaneously 
all the observed diagrams (within the known uncer- 
tainties and biases) both from the point of view of 
the location of the sequences and of the density of 
points in each zone of the diagrams. We start with 
the simplest physically reasonable models and pro- 
ceeded step by step, justifying the increase in com- 
plexity by the requirement of fitting the data. 

5. RESULTS 

Fig. 2 shows that the classical modeling with 
non-evolving homogeneous spheres is incapable of re- 
producing the observed slopes in panels a, b and c. 
But the evolution of the exciting stars is not the only 
factor modifying the spectrum of a giant H n region. 

One of the most striking observed trends is the 
steady increase of [O i] A6300/([O n] A3727+[0 ill] 
A5007) as £W(H/3) decreases. As known, stellar 
winds sweep the surrounding matter and confine the 
emitting gas in a thin expanding shell. This induces 



an important decrease of the ionization parameter 
and allows to reach the highest observed values of 
[O i] A6300/([O n] A3727+[0 in] A5007) in about 
5 Myr (see green curve in Fig. 3). However, it does 
not reproduce the observed slopes in panels a, b, c. 
Also, the models give higher EW(H/3) than observed 
in early ages. Selective dust absorption, if existing, 
would have only a minor effect since the extinction 
is small. The presence of an older, non ionizing stel- 
lar population, attested by stellar features, increases 
the continuum and thus reduces EW(Hf3). However, 
the older population is independent of the age of the 
most recent starburst (although it may depend on 
the metallicity, being likely more important for H n 
galaxies having reached a higher metallicity). The 
blue curve in Figs. 3 and 4 represents a sequence 
with an average old population added. While this 
brings the highest EW(H(3) in better agreement with 
the observations, this does not solve the problem of 
the slopes in panels a, b and c, and does not ex- 
plain why He I A5876/H/3 remains constant in the 
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Fig. 5. The "intermediate" metallicity bin as in Fig. 4. The cyan curve corresponds to the same sequence as the red 
curve in Fig. 4 but with X-rays added. The purple curve corresponds to the same sequence but for a two-component 
chemical composition of the gas. 



entire range of EW(R0). What is needed is to lower 
EW(H{3) more rapidly. One way to achieve this is 
to consider a covering factor of the ionizing source 
decreasing with time. 

While this can explain the "high" metallicity bin, 
the "intermediate" and "low" metallicity bins show 
further problems. The He n A4686 nebular emis- 
sion occurs too frequently and in too wide a range of 
EW(H(3) to be attributed to either the hard photons 
from Wolf-Rayet stars (Schaerer 1998), or the X-rays 
from young starbursts (Cerviho et al. 2002). The 
cyan curve in Fig. 5 shows the effect of photoioniza- 
tion by an X-ray source of luminosity 4 x 10 40 erg s _1 
(10 times the initial H/3 luminosity). Such X-rays 
could originate from massive binaries or from super- 
nova remnants in the age range of 10 - 50 Myr (Van 
Bever & Vanbeveren 2000). Even this amount of X- 
rays cannot heat the gas enough to reconcile models 
and observations in panels d and e. Additional heat- 



ing mechanisms must be investigated. Hard photons 
from planetary nebulae and white dwarfs from earlier 
generations of stars are certainly not enough. Photo- 
electric heating by dust grains is likely negligible, due 
to the low ionization parameters involved. Shocks 
produced by winds and supernovae from the latest 
burst of star formation are far from sufficient. But 
the hypothesis of cither turbulent heating or shock 
heating due to stellar winds from previous genera- 
tions of stars are promising and need to be investi- 
gated. However, a heating deficiency is not the only 
way to explain the discrepancy Points in line ratio 
diagrams can also be displaced by the presence of a 
diffuse gas component. Another option is to invoke 
chemical inhomogeneities, which are indeed expected 
in zones affected by mass loss and supernova explo- 
sions from massive stars. The question is what is the 
state of the newly synthesized matter, and how much 
of it has escaped the nebula. The problem is far from 
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settled. Direct evidence for self-enrichment is scarce 
(see Kobulnicky 1999) but the observed increase of 
[N n] A6584/[0 n] A3727 as £W(H/3) decreases ar- 
gues for a nitrogen self-enrichment on timescales of a 
few Myr. Of course, this in itself docs not argue for 
any oxygen self-enrichment, since the nitrogen and 
oxygen donors are not the same. As an example, 
the purple curve in Fig. 5 combines two sequences 
with nebular metallicities 0.2 Z Q and 0.05 Z Q (the 
apparent metallicity of these combined models as 
derived from temperature based methods would of 
course be intermediate between these values) and is 
good agreement with the observations. Similar con- 
clusions are drawn for the "low" metallicity bin. 

6. FINAL REMARKS 

Considering a large and homogeneous sample of 
H II galaxies with oxygen abundances determined 
from temperature-based methods, we have built 
three subsamples of different metallicities and com- 
pared each of them with models of H n regions ion- 
ized by evolving star clusters of appropriate metal- 
licity. This allowed us to distinguish evolution and 
abundance effects in the emission line sequence of 
H II galaxies. This, and the use of weak lines in 
conjunction with classical strong lines allowed us to 
sharpen up our view on the evolution of H n galax- 
ies. H II galaxies draw a tight sequence not only in 
metallicity but also in age, on a time scale of a few 
Myr, corresponding to the lifetimes of massive stars. 
This is in agreement with the view that the exci- 
tation of an H II galaxy is dominated by a cluster 
of coeval massive stars. We have shown that un- 
veiling the evolution of H n galaxies requires taking 
into account the dynamical evolution of the gas. The 
adiabatic expanding bubble theory provides a good 
frame. While the equivalent width of H/3 is definitely 
linked to the age of the most recent starburst, the 
observations demand that EW(H./3) decreases faster 
than predicted by normal photoionization models. 
This cannot be the effect of an older stellar popula- 
tion alone or of selective dust absorption, and calls 
for a covering factor decreasing with time (implying 
that a growing amount of Lyman continuum radia- 
tion is escaping from H n galaxies). 



The origin of the nebular He n A4686 emission is 
not clear. The question of the heating of giant H II 
regions is not definitely settled. We see evidence for 
an enrichment in nitrogen on a timescale of a few 
Myr. Possibly, H II galaxies are not homogeneous in 
oxygen abundance. 

While our study has allowed to considerably im- 
prove our diagnostic on the evolution of H n galaxies, 
a real understanding of this evolution now implies a 
dynamical modeling that would be able to reproduce 
all the observed trends. 
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